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The non-classical property of subpoissonian photon statistics is extended from one to two-mode 
electromagnetic fields, incorporating the physically motivated property of invariance under passive 
unitary transformations. Applications to squeezed coherent states, squeezed thermal states, and 
superposition of coherent states are given. Dependences of extent of non-classical behaviour on the 
independent squeezing parameters are graphically displayed. 



I. INTRODUCTION 

Non-classical properties and effects of radiation have received considerable attention in the past two decades 
and continue to be an active area of research. Quadrature squeezing [^), subpoissonian photon statistics (SPS), and 
antibunching Q of photons are three prominent examples of such properties leading to measurable effects. Quadrature 
squeezing is related to the reduction of noise in one of the two quadrature components below the coherent state value, 
and has been both theoretically and experimentally studied for one-mode as well as multi-mode fields. Antibunching 
arises when the photon number distribution becomes subpoissonian leading to anticorrelation in the photons detected 
in a typical detection experiment. In all these cases, the diagonal coherent state description of the fields involved does 
not have a classical interpretation and hence no classical description can explain these effects. The extension from 
one to two or more modes for the case of quadrature squeezing is nontrivial and leads to new physical effects ||. 
The phenomenon of SPS has been formulated || and observed j?J primarily for one-mode situations. For situations 
involving two or more modes, in the existing literature, such properties are invariably studied for one of the modes 
or a predefined linear combination of the modes ^ . Such an analysis cannot be used to make any clear statement 
about the classical or non-classical nature of the field involved, because the linear combination of modes which may 
show SPS may in general be different from the mode chosen for the analysis. Another kind of generalisation of SPS to 
two-mode fields has been done using a particular inequality involving the correlation between the two modes; however 
this does not exhaust the possibilities available at the level of quadratic expressions in photon number Q . This clearly 
indicates the need for a more satisfactory way of looking at non-classical statistics, for fields involving two or more 
modes. 

Our aim in this paper is to develop a notion of SPS for two-mode fields which is intrinsically two-mode in character, 
can be used in an unambiguous way to make a statement about the classical or non-classical nature of the field, 
and has physically reasonable invariance properties. The group of linear, homogeneous, canonical transformations 
<Sp(4,3fi), the symmetry group basic to the quantum mechanical description of the two- mode field, naturally splits 
into two parts: the photon number conserving (maximal compact) passive subgroup U(2), and the photon number 
nonconserving (non-compact) active part. The maximal compact subgroup C/(2), while acting on the Hilbert space 
of the two-mode system through its unitary representation, is incapable of generating a non-classical(classical) state 
starting from a classical (non-classical) one because the diagonal coherent state distribution function is covariant 
under such transformations. Therefore, it is reasonable to require that any signature of non-classicality for a two- 
mode system, in particular SPS, be U(2) invariant. To achieve this we regard all modes related to the original ones 
by passive U(2) transformations as basically equivalent; then a survey of the SPS properties for each mode in this 
equivalence class of modes leads to the formulation of a U(2) invariant definition of SPS. We search over the set of all 
modes for that one which minimizes the relevant parameter measuring number fluctuation minus the mean. In this 
way, we arrive at that U (2) combination of the two modes which is most likely to be manifestly subpoissonian. A 
much wider class of non-classical states can be explored using this formalism compared to the earlier ways of handling 
two-mode situations. 

The material in this paper is arranged as follows: In Section |l| we recapitulate the basic kinematics of two-mode 
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systems and the action of the group Sp(A, 3?) on the nonhermitian annihilation and creation operators. The hermitian 
generators of this action, and the maximal compact subgroup U(2) £ Sp(4,$t), are recorded. The notion of U(2)- 
invariant SPS is then developed by regarding all modes related to one another by (passive) 577(2) transformations as 
equivalent, and by minimizing the variable one-mode Q parameter over the group SU(2). The alge brai c machinery 
needed to carry this out, for an arbitrary given state of the two-mode system, is set up. In Section III we consider 
three applications: squeezed coherent states, squeezed thermal states and a general superposition of two coherent 
states. In each case the analytic work is carried out as far as possible, and then we resort to numerical studies which 
are graphically displayed. Section IV contains concluding remarks. 



II. U(2) INVARIANT DEFINITION OF SUBPOISSONIAN PHOTON STATISTICS FOR TWO-MODE 

SYSTEMS 



We consider two orthogonal modes of the radiation field, their orthogonality being achievable by their having 
different frequencies, orthogonal polarizations or different directions of propagation. These modes can be quantum 
mechanically described by photon annihilation operators a r and corresponding photon creation operators a\, where 
r = 1, 2. These operators can be arranged as a column vector £( c ); 



(c) = (d c) ) = 



at 

\4J 



a = 1---4. 



(2.1) 



The superscript (c) on £ indicates that the entries here are complex i.e. nonhermitian. The quadrature components of 
these operators, which are the hermitian phase space variables g's and p's, can be written as another column vector, 
related to £^ c ) by a fixed numerical matrix O: 
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The canonical commutation relations obeyed by the creation and annihilation operators can be written in terms of f; 
or £( c ) : 
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(2.3) 



A general real linear homogeneous transformation on the q's and p's which preserves these commutation relations 
is described by a 4 x 4 real matrix S obeying the condition: 



S p S T = p. 

This is the defining property for the elements of the non-compact group Sp(4, 5i): 



sp(4,n) = { 



= < S = 4 x 4 real matrix 



3 p S T = /3} 



(2.4) 



(2.5) 



When £ undergoes a transformation by S G S'p(4, 5R), the nonhermitian operators £*- c ) transform through a complex 
matrix S^- c \ obtained from S by conjugation with Q: 



S G Sp(A, ft) : C = S( =► 

^(c)/ = S (c) ^( C ) j 



(2.6) 
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The complex matrices are a faithful representation of the real matrix group <Sp(4, 5ft). In this sense we will treat 
them as elements of Sp(4, 5ft). 

The maximal compact subgroup U (2) of Sp(A, 5ft) can be identified as follows: 



£ = 17(2) = {S< c) ([/) S S*p(4,5ft)|£7 e 17(2)} C 5p(4,5ft) 

J7 
U* 

The block diagonal form is responsible for the fact that such transformations do not mix a and a^; in fact JC is the 
largest subgroup with this property. 

Let Tt be the Hilbert space on which £ and act irreducibly. It follows from the Stone- von Neumann theo rem ficfl 
that, since the canonical commutation and hermiticity relations are invariant under the tr ansf ormation (2.6) for any 
€ 5p(4, 5ft), it is possible to construct a unitary operator U{ c -*) on 7i implementing (2.6) via conjugation: 

S {c) e Sp(4,5R) : S$£ c) = U(S^)~ 1 ^U(S^), 

U{S (c) )^U{S {c) ) = 1 on W. (2.8) 

The generators of the operators U(S^) are given by ten independent, hermitian, quadratic expressions in a r and a) r . 
We define the four photon number conserving generators Jo, Jj and the six photon number non-conserving generators 
A,- /•,-./ 1,2,3: 

Jo = \{N + 1) = -(ojoi + a\a 2 + 1); (2.9a) 

Jl = \ ( a l a 2 + 4 a l)> 

J2 = -^i a l a i - a\a 2 ), 



J3 2 



i(aJai-<4o2); (2.9b) 
1 



X a = j(a\ 2 + a 2 1 -4 2 -d$), 
K 2 = -~(4 2 -ai + 4 2 -al), 

K 3 = --(44 + a 1 a 2 ); (2.9c) 

Li = \(a\ 2 -4-4 2 + 4), 

L 2 = \{a\ 2 + a 2 + 4 2 + 4), 

L 3 = --{4 a \ - a i a 2}- (2.9d) 
These generators obey the commutation relations 

[Jj, Jfc] = iejkiJi, 

[J ,J j }=0; (2.10a) 
[Jj,K k or L k ] = ie jk i(Ki or L t ), 

[J , ± iL 3 ] = T(Kj ± iij); (2.10b) 
[isT,-,2ffc] = [7j,7 fe ] = -iejkiJi, 

[Kj,L k ]=iS jk J . (2.10c) 

From the above commutation relations, it is clear that Jo and Jj form the algebra of U (2) and hence generate the 
unitary operators corresponding to the elements of the maximal compact subgroup K, of Sp(A, 5ft). On the other hand, 
Kj and Lj are the generators of the unitary operators corresponding to the non-compact elements of Sp(A, 5ft) and 
they do not form a closed algebra. These non-compact elements are actually the squeezing transformations and their 
complete classification has been given elsewhere [||. 
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We now consider the notion of SPS for the physical states of a two-mode system. For one-mode systems, such an 
analysis is simple and is based on Mandel's Q parameter Q: 

q= W -<•'■>' ,2.1, 
{a'a) 

where a and are the annihilation and creation operators for the one-mode radiation field, the expectation values 
being taken for the state of interest. The Q parameter distinguishes between physical states as having poissonian, 
subpoissonian and superpoissonian photon statistics, as Q is 0, < and > for the above cases respectively. In 
particular, the states with negative Q are non-classical, in the sense that such a distribution can not be derived from any 
classical statistical ensemble. Therefore, in this limited sense, the Q parameter can be used to classify states as classical 
and non-classical. More precisely, Q < (> 0) is a sufficient (necessary) condition for non-classicality(classicality). 

For a situation involving two modes, the notion of SPS defined above is not appropriate. At the most, one can 
analyze the photon statistics of one of the modes, or a preselected linear combination of both. Then again, for a given 
state, this mode which one chooses need not be the one in which the photon number distribution may be non-classical. 
Hence the sign of Q for the preselected mode may not disclose the non-classical nature of the two-mode state, even if 
it is non-classical. This clearly indicates that an intrinsically two-mode notion of SPS is required. 

The standard way O of distinguishing the classical from the non-classical states (already implicitly assumed in 
the above) is through the diagonal coherent state description. The general two-mode coherent state with complex 
two-component displacement z — (2:1,2:2) is defined by 

\z) = exp (z ■ a) - S* ■ a) |0, 0) 

= exp^-i|zi| 2 -i|z 2 | 2 ^ exp (z ia \ + z 2 <4) l°>°>- ( 2 - 12 ) 

These are normalized states and form an over-complete set. A given two- mode density operator p can be expanded 
in terms of them: 

Zl Z2 <p(zi,z 2 )\zi,z 2 ){zi,z 2 \ (2-13) 

7T" 2 

The unique normalized weight function <p(zi, z 2 ) gives the complete description of the two-mode state and can in 
general be a distribution which is quite singular fll|| . In the case when ^(2:1,2:2) can be interprete d as a probability 



distribution (i.e. it is nonnegative and is nowhere more singular than a delta function), equation (2.13) implies that 
the state p is a classical mixture of coherent states which have a natural classical limit. Such quantum states are 
referred to as classical; in contrast the others for which (f>(zi,z 2 ) either becomes negative or more singular than a 
delta function somewhere, are defined as being non-classical. This classification is general and can be done for any 
number of modes. In particular, for the one-mode case, the states having negative Q are a subset of the states with 
non-classical diagonal coherent state distribution functions. 

When the two-mode state, with density matrix p, transforms under a unitary operator corresponding to the compact 
U(2) subgroup of Sp(4, 5i), the distribution </>(zi,Z2) undergoes a point transformation given in terms of the 17(2) 
matrix U E U{2): 

p' = uis^mpU^Xu))- 1 O 0'(*i,«O = 4>(z[,z' 2 ), 

= U^) (2.14) 

Thus, under U(2) transformations classical states map on to classical ones and non-classical states to non-classical 
ones; these transformations are incapable of generating a non-classical state from a classical one. Therefore, it is 
reasonable to demand that any signature of non-classicality be invariant under such transformations. 

At this stage, we recapitulate and collect some interesting and important properties of the maximal compact 
subgroup K. of S'p(4,3?): 



(a) As is clear from eqn.(2.7), when £( c ) undergoes a U(2) transformation, the annihilation operators a r 's are not 
mixed with the creation operators aps. 

(b) The action of the elements of U (2) (generated by Jo and Jj) on a state does not change the total photon number 
or its distribution. 
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(c) The diagonal coherent state distribution function is covariant under U(2) transformations. 

(d) One requires only passive optical elements to experimentally implement any U(2) transformation on a state of 
the two-mode electromagnetic field |l2[ . 

Motivated by the above considerations we now define an intrinsically two-mode and U(2) invariant notion of SPS. 
For the purpose of our present analysis it is convenient to define U(2) transformed mode operators in terms of two 
column vectors A and a: 

*=(:t), «=(*) (2.i5) 



° 2 / V a2 

where a\ and a 2 are complex numbers such that: 

U W = (-L ai) eSU{2) > M 2 + M 2 = 1 ; 

U(a, if,)=( _^J oa Jft^ \ e U(2), < V < 2tt (2.16) 

When £( c > undergoes a U(2) transformation given by U(a,ip), the annihilation and creation operators for the first 
transformed mode can be written in terms of A and a alone: 

a(a) = a^A — a\a\ + 0^0-2 

a(a)^ = A^a = a.\a[ + a^aJ,. (2-17) 

Thus the most general normalized "first mode" after the U(2) transformation is determined by SU(2) € U(2) inde- 
pendent of ip. This particular mode will henceforth be called the SU{2) transformed mode, and a will be used to 
denote the SU{2) element involved. 

Let p be the density matrix for any (pure or mixed) state of the two-mode radiation field. Then we can define the 
following function: 



a(a) t2 : a(a) 2 ) p - (a(a) t a(a))p 

(2.18) 



Tr{pa{ a y 2 a{a) 2 ) - (Tr(pa(a)U(a))) 2 



TripA^A) 

which is similar to the Mandel Q parameter for the SU{2) transformed mode a{a). 

When the state p is transformed by the unitary operator U{S^ C \U)) for some U € U(2), the function Q(p; a) can 
be shown to change covariantly: 

S (c) (t/)G/C : p' =U{S {c) {U))pU(S {c) (U))- 1 => 

Q(p'; a) = Q{p; a'), a' = Ua (2.19) 

Now an overall phase change corresponding to elements in the U(l) subgroup of U{2) actually leaves Q(p;a) un- 
changed, therefore no dependence on ip has been shown. So we have the freedom of running over all a's € SU (2) i.e. 
we can choose various linear combinations of the two modes involved, related to each other by SU (2) transformations. 
Since we want to unearth the signature of the non-classical nature (if present) of the photon statistics, we vary a till 
we reach the minimum value of the function Q(p; a): 

Q{p)= Min Q{p;a) = Q(p;a) s.t. Q(p;a)<Q(p;a) (2.20) 
over all 

a£SU(2) 

If Q(p) < we conclude that the photon number distribution for the two- mode state p has a non-classical feature and 
we call it subpoissonian, or amplitude squeezed. This is our U(2) invariant definition of SPS for states of two-mode 
fields. The mode in which the subpoissonian nature is manifest to the maximum degree is a(a). 

The numerator in our definition of Q(p\ a) consists of two terms, one arising from the expectation values of quadratic 
expressions in the creation and annihilation operators and the other arising from the expectation values of quartic 
terms. The quadratic term can be written: 
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Tr (pa (ay a (a)) = s + u.q, 

q = q(a)=a t &a, (2.21) 

with the dependence on the state p and on a £ SU (2) being clearly separated. The state dependent variables s and 
u transform under SU(2) like a scalar and a cartesian vector respectively, and can be evaluated from the equation 

Tr{pa\.a s ) = s8 rs + Uj{(Tj) rs r,s = 1,2. (2.22) 

The term involving the expectation values of quartics in a r and a\ can be written in terms of the non-compact 
generators K and L of Sp(A, 5ft), and a vector A representing the SU(2) element involved: 

2 1 

Tr(pa(ay a(a) 2 ) = -AjA£i?, fe , 

H jk = H* kj = TripiKj - iL^Kk + %L h % j, k = 1,2, 3 , 
A = A (a) = -ia T a 2 va, A (a). A (a) = 0. (2.23) 

The hermitian matrix can be written in terms of two real matrices, the real symmetric R and real antisymmetric 
S, as _ff = R + iS. The matrix R transforms under SU(2) as a second rank tensor whereas the matrix S can be 
represented by a cartesian vector v under SU(2), related to S by Vj = hejklSkl- 

The denominator of Q(p; a) is [7(2) invariant since the operator „4 T ./4 = a|ai + a\a2 is t/(2) invariant; it docs not 
depend upon a and can be written in terms of s as: 

V(Q(p; a)) = Tr(pA^A) = 2s (2.24) 

After some algebra, the complex vector A can be eliminated in favour of the real vector q, and Q(p; a) can be written 
in terms of the state dependent symmetric second rank tensor _R, the vectors u, v and the scalar s as: 

Q(p; a) = Q{p; q(a)) = j^(TrR - q q k R ok + 2i.q - 4(s + u.q) 2 ) (2.25) 

Using the U(2) covariance of a), we can assume without loss of generality that the real symmetric matrix R is 
diagonal, and eq( [2.25| ) then takes the simpler form: 

Q(p;q(a)) = ^-(TrR-Y,^Rj 3 +2v.q~A(s + u.q) 2 ) (2.26) 

3 

The dependence of Q(p;q(a)) on a € SU(2) is through the real unit vector q(a), which can be represented on the 
surface of a unit sphere. In order to obtain the invariant Mandel parameter Q(p) for a given two-mode state, we have 
to minimize Q(p; q(a)) with respect to q(a), the parameters R, v, s,u being determined by p. The most convenient 
coordinates which one chooses on the surface of the sphere to carry out this minimization will depend upon the 
physical state p under consideration. 



III. APPLICATION TO TWO-MODE SQUEEZED COHERENT STATES, SQUEEZED THERMAL 
STATES AND SUPERPOSITION OF COHERENT STATES 

In this Section, we apply the formalism developed in the previous Section to various interesting two-mode states. 
Here we will see the relation with the classification of two-mode squeezing transformations given in . 



A. The case of squeezed coherent states 



The most general squeezed coherent state is obtained by applying the operator U(k,l) = e l ( k - K + LL ) to the two- 
mode coherent state 1 21,2:2) defined in cqn. ( 2.12| ), for some complex 21,22, where K and L are the non-compact 
generators of Sp(4, K) defined in eqn. (2.9) and k and I are real vectors. The operator U(k,l) is conjugate to 
V(( '(a, b) = expi(aK 2 + bL\) for some a > b > 0, via an operator U(S C (U)): 
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U(k, I) = U- 1 {S (c) {U))U {0) (a, b)U(S {c) (IT)), 
U^b) = exp (<iZ*> (af - a?)), exp (<i+*> (of - a|)) (3.1) 



Each U^\a,b) is a representative of an equivalence class of two-mode squeezing transformations. For a = b we have 
the essentially single mode case, while for b = we have maximal involvement of the two modes. For the minimization 
of 1/(2) covariant Q(p; a), the overall 17(2) factor U~ x (S^ C \U) is irrelevant. Also, the action of the operator U(S^ (U)) 
on ^1,22) transforms it into another coherent state \z[,z 2 ), with z[, z' 2 related to 21,2:2 through the corresponding 
17(2) transformation. Thus it suffices to examine the particular class of squeezed coherent states 

\ Zl ,z 2 ,a,b) =UW(a,b) \z u z 2 ). (3.2) 

A complete discussion of the two-mode squeezing transformations and squeezed states has been given in [Q . 

The Mandel parameter Q(zi, z 2 ,a,b;q(a)) for the SU{2) transformed mode for squeezed coherent states can be 
calculated by str aigh tforward algebra and turns out to be rather lengthy. The complete expression is given in the 
appendix (eqn. ( |A1[ )). Q(z\, z%, a, b; q(a)) depends on a, b through hyperbolic functions and on |zj.|,|z2| through 
polynomial functions. Its dependence on the phases of 21 and 22 and the polar coordinates on the surface of the 
unit sphere describing the unit vector q(ot), is through trigonometric functions and is oscillatory in nature. In order 
to obtain the invariant Mandel parameter, this function has to be minimized with respect to q(a). Since this is not 
possible analytically, the results obtained numerically are displayed in Figures ([!]) , ^) and (||) jl3| . 

In each figure, we plot the minimum value of Q(zi, z 2 , a, b; q{a)) as a function of the squeeze factors a and b, keeping 
the complex displacements 21 and 22 fixed. Figure ([l](a)) displays the results for the squeezed vacuum; this never 
shows SPS. The plots of Figures (|l](b),(c) and (d)) on the other hand are obtained by varying the phase of one of the 
displacements (22), keeping its magnitude fixed, with the other displacement (21) being zero. Different values for the 
phase of the non-zero displacement give qualitatively different results; in particular when this phase is |_, as is clear 
from Figure (0(d)) even some of the essentially single mode states lying along a — b show SPS. In Figure (g) we choose 
equal magnitudes of displacements for the two modes; plots have been generated for different values of their phases. 
The displacement parameters in Figure @ are unequal in magnitude; four plots have been given for the same choices 
of phase values as in the corresponding plots in Figure (|^) . The qualitative features of individual plots are similar to 
the corresponding plots in Figure (||) though the actual values of the invariant Mandel parameter are different. 

We now make some general remarks about the results described above. In all the plots of the three Figures (Q), 
(|J) and (||), every point in the region b > a can be mapped onto a corresponding unique point in some region a > b 
(which is in general not in the same figure), through that 1/(2) transformation of the displacements z\ and z 2 , which 
effectively changes U^(a, b) to U^(b : a). Whenever the displacement parameters are invariant under this particular 
1/(2) transformation, the plot has a symmetry about the line a = b; as in all the plots of Figure (0). Such a symmetry 
is not exhibited by the plots of Figures (|^) and (|J) . In all the plots the invariant Mandel parameter is zero or negative 
along the line a = b i.e. for the subset of essentially single mode squeezed states. This happens because, even though 
the choice of displacement parameters is such that the single mode which is squeezed has superpoissonian statistics, 
(Q > 0), the minimisation chooses the other mode which is in a coherent state (Q = 0). Apart from the case of 
squeezed vacuum(Figure (|l|a)) all other choices of displacement show SPS for some values of the squeeze parameters 
a, b. When squeezing becomes large in comparison to the displacement, and we are away from the line a = b, SPS 
disappears and the states tend to become more and more superpoissonian. 

B. The case of squeezed thermal states 

We next look at the case of a two-mode isotropic thermal state subjected to squeezing. The normalized density 
operator corresponding to the inverse temperature f3 = huj/kT is explicitly C/(2) invariant and described in the Fock 
representation by: 



Po 



(/3) = (1 - e- 13 ) 2 exp [-(3(a\ a a + a\ a 2 ) 



= (l-e-P) 2 J2 e^ (ni+n2 Vi»«2)<ni,ri 2 |, (3.3) 

ni,n 2 =0 



with U (2) invariance expressed by: 



e idJ ° potf) e~ m ° =e ia ' J p o (0) e~ ld ' J = p {f3). (3.4) 
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Therefore it suffices to examine the properties of the density operator obtained by conjugating po(3) with l(( \a, b): 



p(3; a, b) = «(»»(«, b) p {8) W(°)(o, ft)" 1 . (3.5) 
The Mandel parameter Q(/3; a, b; q(a)) for the state p(8; a, b) is calculable by straightforward algebra: 
Q{8;a,b;q(a)) = 

[{e - 1) (2 (1 - e p ) + 2 (1 + e p ) cosh(2a) cosh(26))] _1 x 

~ ((1 - q 2 ) (2 (1 - ef 3 ) 2 +4(1- e 2 ?) cosh(2a) cosh(26) + (1 + e* 3 ) 2 ( cosh(4a) + cosh(46))) 
+ (1 + e' 3 ) 2 ( 9 i 2 - g 2 2 )(cosh(4a) - cosh(46)) 

~\ (1 + 93 2 ) (10 ~ 12e' 3 + 10e 2/3 + 16 (1 - e 2/3 ) cosh(2a) cosh(2&) + 6 (1 + e 13 ) 2 cosh(4a) cosh(46))^ 

--((1 + e 13 ) q 3 (4 - Ae 13 + 6(1 + e 13 ) cosh(2a) cosh(2fe)) sinh(2a) sinh(2fe)) 
1 



-((2 - 2eJ 3 + 2(-l + e 2p ) cosh(2a) cosh(26)) + ((1 + e 13 ) q 3 sinh(2a) sinh(26))) 2 



(3.6) 



here qi,q 2 , 93 are the cartesian components of q with q\ + q\ + q\ = 1. 

The minimum value of the function Q{j3,a,b,q{a)), the parameter Q(p(/3, a,b)), can be calculated analytically. 
The state p(/3, a, b) being the squeezed thermal state is always superpoissonian. For a given temperature (given (3) 
this superpoissonian nature is least for the case when only one mode is squeezed(a = b), increases as the squeezing 
becomes increasingly two mode in nature and finally is maximum when the state is maximally two mode squeezed i.e. 
when a = (b = 0) for a given b (a). When the temperature is changed the states with higher temperatures (lower 0) 
are more superpoissonian compared to the ones at lower temperatures (higher /3).Thus for fixed a and b, Q(f3;a,b) 
increases as 3 decreases. The actual plots of Q(p(B,a,b)) as a function of a, b are given at different temperatures in 
Figure(|) II . 

It is interesting to note that the particular mode for which the function Q(p(3,a,b),q(a)) is minimum turns out 
to be one of the original modes, corresponding to q 3 = ±1. This happens because the thermal state density matrix 
Po(3) is explicitly U(2) invariant and the representative two-mode squeezing op erat or (a, b) can be factorized into 
two commuting operators, each pertaining to one of the original modes see eq. ([O]). In general, for a different choice 
of the representative operators, the minima could occur at an arbitrary SU (2) transformed first mode. All the plots 
of Figure( ||) are sy mm etric about the line a = b because of the explicit U (2) invariance of the thermal state density 
matrix po(p) (eqn.( 3.4)). 



C. The case of superposition of coherent states 



Lastly we apply our formalism to the superposition of two two-mode coherent states. In this case, no squeezing 
transformation U^(a, b) is involved. For simplicity we consider only the case with real displacements. 

A general superposition of two two-mode coherent states with real displacements and a phase difference r\ between 
them is given by: 

\ip(u 1 ,u 2 ,v 1 ,v 2 ,r,r])) = — (|ui,u 2 ) + r cxp (irj)\v 1 ,v 2 )) 

where iV 2 = 1 + r 2 + 2 r cos r\ exp ^— i(u 2 + u 2 + v\ + v 2 ) + u 1 v 1 + u 2 v 2 ^j (3.7) 

With the help of a U(2) transformation, without any loss of generality we can set v 2 — and thus it suffices to 
study only the states \tp(ui, U2, v 1, 0, r, 77)). The Mandel parameter for the SU(2) transformed mode, the function 
Q(ui, u 2l Vi, r, rj\ q(a)) for this superposition of coherent states is given in terms of the polar coordinates 8 and <f> on 
the surface of the sphere representing q as: 

Q(ui,U2,Vi,r,rj;q(a)) = 
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4 (uf + u\ + r 2 v\ + 2 e -^«i 2 +«2 2 +«i 2 )+«i v x r Ui Vi CQS ^ N 
(l + r 2 + 2 e -^(^ 2 +(«i-fi) 2 ) r cos(77)) x 

/ q 4 ^ 4 ^ 2 

( 4 (ui 4 + r 2 vi 4 ) cos(-) +4u 2 4 sin(-) + 8 ui 3 u 2 cos(-) cos(</>) sin(6») 



+8uiu 2 3 cos(</>)sin(-) sm(9) + 2 u 2 u 2 (2 + cos(2 0)) sin(0)' 
+r (V^2 2 +(«i-«i) 2 ) (l + r 2 ) +2e- U22 - ( - U1 -" 1 ') 2 r cosfa)) x 



8 ui 2 «i 2 cos(-) cos(r7) 
2 

+8 ui u 2 «i 2 cos(-) cos(?7 + </>) sin(0) + 2 u 2 2 t>i 2 cos(?7 + 2 <^>) sin(6 1 )' 

/ Q 2 ^2 

f 2 (ui 2 + r 2 «i 2 ) cos(-) +2u 2 2 sin(-) + 2 ui u 2 cos(0) sin(6») 

1 1 2 2 / 9 2 

+ _ e -^("2 2 +2 («i-«i) 2 ) r / 4 Ui Wi cos (_) cos(?7) + 2 u 2 ui cos(?7 + 4>) sm{8) 



(3.8) 



The minimum values of this function with respect to 9 and (j) have been computed numerically and the results are 
shown in Figure (|). Each plot in this figure contains two curves showing Q{p) as a function of the relative phase r\ 
corresponding to two different values of relative weight factor r. The amount of SPS varies with the relative phase 
in a similar way for all the plots. For all parameter values in all plots Q(p) < 0. This happens because the most 
general superposition of two two-mode coherent states can be transformed with the help of a U (2) transformation 
into a product state with one factor being a coherent state, and the other a superposition of two one-mode coherent 
states: 



i(K) + r exp («7)K»K) = ^U{S^{U)){\u u u 2 ) + r cxp (z ?7 )K^» 



(3.9) 



Thus when Q(p; a) turns out to be nowhere negative the minimization chooses that U(2) transformed mode which is 
in a coherent state. 

It is i nter esting to point out that for a factorized two-mode state such as the expression on the left hand side 
of eqn. (|3.9| ), the mode a(a)which minimizes Q(p;a) is generally neither of the two initial modes but a nontrivial 
combination of them. 



IV. CONCLUDING REMARKS 



The main aim of this paper has been to develop a specific signature of non-classicality for two-mode states. Both 
quadrature squeezing and SPS are well defined concepts for a single mode. In this paper we have extended the notion 
of SPS to two modes by showing how to choose the appropriate single mode which shows SPS to the maximum extent, 
considering all modes related to each other by passive U(2) transformations as equivalent. A similar treatment of 
quadrature squeezing has been given elsewhere. 

We would like to emphasize the subtle role played by the choice of the denominator of Q(p; a). Any choice which 
is everywhere non-negative will not change the qualitative results obtained from the minimization of Q(p; a) i.e. the 
super or subpoissonian nature of the state p. However the extent of SPS, and the location of the most non-classical 
mode, depend upon the exact choice one makes for the denominator. To illustrate this point we choose the two- mode 
Fock state \ni,n 2 ). The Mandel parameter for the SU{2) transformed mode is given by: 

Q(ni,n 2 ,q(a)) = — — x (-2 (m + n 2 ) + (m + n 2 ) 2 + 

4 (m + n 2 ) 

(ni(l - ?ii) + «a(l - n 2 ) ) ( qi 2 + q 2 2 ) - 2 (m - n 2 ) q 3 - (m + n 2 f q 3 2 ) (4.1) 
This function reaches its minimum at 03 = +1 with minimum value — „ j 1 ^ for n\ > n 2 and at q% = — 1 with the 

I L\ T li 2 

minimum value — — ^ for n 2 > n\. Thus for our U{2) invariant choice of the denominator Tr(pA^A), for a Fock 
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state, the mode with the larger number of photons is more non-classical. On the other hand if one chooses the U(2) 
covariant denominator Tr(pa(ay a(a)) , for a Fock state, both the modes are equally non-classical irrespective of the 
number of photons present in each mode: the minimum value of this alternatively defined parameter is —1 for each 
mode. 

For one- mode fields the Mandel parameter can be written as a function of the number operator a^a and hence 
is determined by (the moments of) the photon number distribution. In contrast, for two-mode fields the Mandel 
parameter for the SU(2) transformed mode can not be expressed as a function of the number operators a\a\ and 
a 2 d2 and therefore is not determined by the photon number distributions in the original modes. There could be other 
signatures of non-classicality which are meaningful at the one-mode level and can be extended in the spirit of this 
paper to more than one mode. In contrast, it will be interesting to explore the possibility of having signatures of 
non-classicality which are not definable at the one-mode level at all, but are present only at the two-mode level. These 
will be presented elsewhere. 
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APPENDIX A: 



We give here the function Q(z\, z 2l a, b; a) for the squeezed coherent state with z\ = ue %Vu and z 2 = ve liPv . The first 
term is the denominator, followed by the numerator terms arranged according to their dependence on a and b. First 
the terms independent of a, b appear, followed by the ones depending upon a or 6 alone, and then the ones depending 
on both a and b. The last three terms originate from quadratic expressions of creation and annihilation operators and 
are not arranged. 

Q(a; zi, z 2 ,a, b) = 

2 [-2 + cosh(2 (a - &)) + 2u 2 cosh(2 (a - &)) + cosh(2 (a + &)) + 2v 2 cosh(2 (a + &)) 
+ 2 u 2 cos(2 <p u ) Sinh(2 (a - b)) + 2 u 2 cos(2 p„) Sinh(2 (a + &))] _1 x 

1(5 + m 4 + m 4 + 2 (u 2 + v 2 )) 

o 

+ 1 (v 4 cos(4p„) (-l + cos(0))) 

8 

+ i(( u 2 - w 2 ) (2 + m 2 + m 2 ) cos(0)) 
8 

-^(u 4 cos(4^„) (l + cos(0))) 



16 

+ 



' '\ + u 2 (2 + u 2 - u 2 cos(4 (p u )) +v 2 (2 + v 2 - v 2 cos(4 p„))) sin(6') 1 



I („,, cn*h(<?n\ ( ~ 8 cos ^ cos^-^) - (2 + u 2 + v 2 + (u 2 - v 2 ) cos(0)) sin(0) sin(p u -p„) \\ m 
4 ^ C0S1 V a J ^ _ u 2 (1 + cos(e)) sin( ^ sin(3 ipu + ^ v)+v 2 ^_ cos(e)) sin( ^ sin( ^ + 3 Vv) J J sin^; 

+ 1 f u „ ( ~ 8 cos (^) cos (^ + <Pv) - (u 2 -v 2 + (2 + u 2 + v 2 ) cos(6>)) sin(» sin(^„ + tp v ) \ Sinh(2a) 
4 y ^ — m 2 (1 — cos(#)) sin(0) sin(y>„ — 3 tp v ) — u 2 (1 + cos(0)) sin(0) sin(3 <p„ — tp v ) J ^ ' ' 

If , , / 8 sin(0) sin(<p u - sin(0) + cos(0) cos((/j„ - ip v ) (2 + u 2 + v 2 + (u 2 - v 2 ) cos(0)) 

4 \ \ - cos(0) 2 m 2 cos(3 yj„ + <p„) cos(|) - 2 w 2 cos(^) cos(<£ u + 3 ip v ) sin(|) sin(0) 

-u 1 f,,,, f 1,2 cos ^ cos(<p u -3<p„) (-1 + cos(0)) +u 2 cos(0) cos(3^ u -<p„) (1 + cos(0)) \ 

+ 4 r W ^ - cos(0) cos(^ + p„) (m 2 - « 2 + (2 + u 2 + m 2 ) cos(0)) - 8 sin(0) sin(^„ + p„) J SmW bmh ^ ^ 

1 6 4 

+ - (3 + 12 u 2 + 6 w 4 + 2 u 4 cos(4 p„)) cos(-) cosh(4(a - 6)) 
8 2 

1 4 

+ - (3 + 12 v 2 + 6 v 4 + 2 v 4 cos(4 tp v )) cosh(4(a + b)) sin(-) 
8 2 



1 / 4 \ 

+ 2 ( u2 (3 + 2m 2 ) cos(2</?„)cos(-) Sinh(4(a - &))J 

1 / 4 \ 

+ - (v 2 (3 + 2w 2 ) cos(2^„)sin(-) Sinh(4(a + &))J 

(^(l + 2u 2 ) cos(^) ^ + 
+ cosh(2(a- &)) |-((1 + 2m 2 ) sin(^ 



., 2 ^ u 2 w 2 cos(2^„) sin(2 0) sin(2^)sin(0) z 

+ cosn(2(a — oj) 1 1,1 



. 6 2 \ u 2 v 2 cos(2^) sin(2 0) sin(2 y„) sin(6») : ' 
2 cos(2 (1 + cos(0))) + V " ' — ; - ^"^^^^ I Sinh(2(a - 6)) 



(1 + 2 ?i 2 ) w 2 sin(2 0) sin(2 <p„) sin(6>) 



2 N 



u 2 (l + 2 w 2 ) sin(2 sin(2 <^„) sin(6»)" 



+ ^w 2 cos(2 (-1 + cos(6»)) i '- '—^ —J Sinh(2(a - &)) 

+ i (^((l + 2u 2 ) (l + 2t; 2 ) +2ii 2 i> 2 cos(2 0) cos(2^ u ) cos(2^)) cosh(2(a-6)) cosh(2(a - &)) sin(6») 2 ^) 
+ i ((4m 2 (1 + 2m 2 ) cos(2<^„) + 2 (l + 2u 2 ) v 2 cos(2 0) cos(2^)) cosh(2(a - b)) sin(9) 2 Sinh(2(a - 6))) 
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+ i(((l + 2?i 2 ) (l + 2v 2 ) cos(2</>) +8u 2 v 2 cos(2 tp u ) cos(2 sin(6>) 2 Sinh(2(a - 6)) Sinh(2(a - 6))) 



+ - ^(2ii 2 (l + 2w 2 ) cos(2 0) cos(2 </?„) + 4 (l + 2u 2 ) w 2 cos(2<p„)) cosh(2(a - b)) sin(0) Sinh(2(a - b)) 

if 9 2 \ 

+ — ( u v cos(0) (3 (l + u 2 ) cos((p u — ip v ) + u 2 cos(3 tp u + ip v )) cos(-) cosh(2(2a — b)) sm(9) J 

if 9 2 

+ — ( u v cos(0) (u 2 cos(3 ip u — tp v ) + 3 (l + u 2 ) cos(</?„ + ip v )) cos(-) sin(#) Sinh(2(2a — b)) 

If 9 2 

+ — iuv cos(0) (3 (l + v 2 ) cos((^„ — tp v ) + v 2 cos {ip u + 3 ip v )) cosh(2(2a + b)) sin(-) sin(0) 

if 9 2 N 

+ - (uv cos(0) (w 2 cos(</?„ — 3 + 3 (l + v 2 ) cos(<^„ + ip v )) sin(-) sin(#) Sinh(2(2a + &)) 

1 / 2 \ 

+ — ( uv cos(-) sin(</>) (— (u 2 sin(3<^„ — + 3 (l + u 2 ) sin(<p u + sin(#) Sinh(2(a — 2 6)) j 

If 9 2 \ 

+ - ( mi sin(0) (w 2 sin(</>„ - 3ip v ) + 3 (l + w 2 ) sin(^„ + <p„)) sin(-) sin(0) Sinh(-2(a + 26)) J 

If 9 2 

— — [uv cosh(— 2(a + 26)) sin(<?!>) (3 (l + v 2 ) sin(< / 9„ — + w 2 sin(y>„ + 3 tp v )) sin(-) sin(0) 

If 9 2 

+ — ( uvcos(-) cosh(2(a — 2 6)) sin(</>) (—3 (l + u 2 ) sin(ip u — ip v ) + u 2 sin(3 (p u + fv)) sin(0) 



+ 



i (u 2 w 2 cos(2</>) sin(2<^„) sin(2 93^) sin(0)' 



+ 



1 (l + cos(0)) ((l + 2u 2 ) cosh(2(a-6)) + 2w 2 cos(2<^)Sinh(2(a-6))) 

2 + 4 

(l-cos(0)) ((l + 2v 2 ) cosh(2(a-6)) + 2v 2 cos(2 ip v ) Sinh(2(a - 6))) 



sin(0) ( Sin( -^' ) ( cosh (^ 26 ) sin (^« _ V«) _ sin (^ + Vv) Sinh(-26)) \ 
^ ' \ +cos(0) (cos(<^„ — ip v ) cosh(2a) + cos(tp u + ip v ) Sinh(2a)) J 
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FIG. 1. Plots of the invariant Mandel parameter Q(p) for squeezed coherent states as a function of squeeze parameters a and 
b. Fig. 1(a) shows the plot for squeezed vacuum i.e. z\ = Z2 = 0. Figures(l(b), (c) and (d)) show the plots for \z\ [ = 0, = 3.0 
and the phase of Z2 taking the values 0, 7r/4 and 7r/2 respectively. 



FIG. 2. Plots of the invariant Mandel parameter Q(p) for squeezed coherent states as a function of squeeze parameters a 
and b for the case when the magnitudes of the displacements in the two modes are equal: [zi| = \z%\ = 2.0. The values of the 
phases of z\ and Z2 in Figures(2(a),(b),(c) and (d)) are (0,0), (0, 7r/4), (0, 7r/2) and (tt/2,tt/2) respectively. 



FIG. 3. Plots of the invariant Mandel parameter Q(p) for squeezed coherent states as a function of squeeze parameters a 
and b for the case when the magnitudes of the displacements in the two modes are unequal: \zi\ = 2.0 and \z%\ = 4.0. The 
values of the phases of z\ and Z2 in Figures(3(a),(b),(c) and (d)) are (0,0), (0, 7r/4), (0, 7r/2) and (n/2,n/2) respectively. 

FIG. 4. Plots of the invariant Mandel parameter Q(p) for squeezed thermal states as a function of squeeze parameters a and 
b at different inverse temperatures; /3 takes the values 0.5, 1.0, 2.0 and 4.0 in Figures(4(a),(b),(c) and (d)) respectively. 



FIG. 5. Plots of invariant Mandel parameter Q(p) for superposition of two two-mode coherent states as a function of the 
phase difference between the two states are shown for two different values of relative weight, r = 0.5 and r = 1.0, and a given 
set of displacements. Values of displacements Ui, U2 and V2 are (0.5,0.5, 1.0), (0.5, 1.0, 1.0), (1-5, 1.0, 1.0) and (1.5, 1.0,0.5) 
for Figures(5(a),(b),(c), and (d)) respectively. 
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